arXiv:1501.07064v2 [astro-ph.CO] 25 May 2015 


A NEW DYNAMICAL MASS MEASUREMENT FOR THE 
VIRGO CLUSTER USING THE RADIAL VELOCITY 
PROFILE OF THE FILAMENT GALAXIES 
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ABSTRACT 


The radial velocities of the galaxies in the vicinity of a cluster shows devi¬ 
ation from the pure Hubble flow due to their gravitational interaction with the 
cluster. According to a recent study of Falco et al. (2014) based on a high- 
resolution N-body simulation, the radial velocity profile of the galaxies located 
at distances larger than three times the virial radius of a neighbor cluster can be 
well approximated by a universal formula and could be reconstructed from direct 
observables provided that the galaxies are distributed along one dimensional fil¬ 
ament. They suggested an algorithm for the estimation of the dynamic mass of 
a cluster M v by fitting the universal formula from the simulation to the recon¬ 
structed radial velocity profile of the filament galaxies around the cluster from 
observations. We apply the algorithm to two narrow straight filaments (called 
the Filament A and B) that were identified recently by Kim et al. (2015) in the 
vicinity of the Virgo cluster from the NASA-Sloan-Atlas catalog. The dynamical 
mass of the Virgo cluster is estimated to be M v = (0.84^g;^) x 10 15 /i - 1 M 0 and 
M v = (3.241^3®) x 10 15 h~ l M Q for the cases of the Filament A and B, respectively. 
We d i scuss observ ational and theoretical systematics intrinsic to the method of 
Falco ef al. 120141) as well as the physical implication of the final results. 


Subject headings: galaxies : clusters — large-scale structure of universe 


1. INTRODUCTION 

The simplest and probably the most reasonable explanation for the current accelera¬ 
tion of the universe is the anti-gravitational effect of the cosmological constant (A) that is 
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supposed to be dominant at the present epoch. The fundamental assumption that under¬ 
lies this explanation is that the general relativity (GR) is the universal law of gravity, the 


validity of which has been confirmed by recent cosmological observations flReves et al.l 12010 


Woitak et al. 2011 : Rapetti et al. 12013 ). Besides, the latest news from the Planck CMB (Cos- _ 

mic Microwave Background) mission has drawn the following bottom line flPlanck 2015 results XIV 
20151 ): The risk that one has to take by accepting the ACDM (A+cold dark matter) cosmol¬ 
ogy founded upon GR is lower than that by accepting any other alternative. 

Nevertheless, the conceptual difficulty to acquiesce in the extremely fine-tuned initial 
conditions of the universe associated with A still dares the cosmologists to take a high risk of 
suggesting modified gravity (MG) models and developing independent tests of GR. One of 


Schmidt et al. 

2009; 

Schmidt 

2010: 

Zhao et al. 

2011) 


cluster with its dynamic mass (e.g., see 
The former is estimated from the weak gravitational lensing effect generated by the mass of 
a cluster on the shapes of the foreground galaxies, while the latter is conventionally derived 
from the velocity dispersions of the central satellites of a cluster. A substantial discrepancy 
between the two mass measurements for a cluster, if found and accurately measured, would 
challenge the key tenet of GR. 

However, it is a formidable task to measure the dynamical mass of galaxy clusters with 
high accuracy. The success of the conventional method of using the satellite velocity dis¬ 
persions is contingent upon the completion of relaxation in the dynamical state and the 
spherical symmetry in the shape. There is another method for the dynamical mass measure¬ 
ment of the clusters that uses the peculiar velocities of the neighbor galaxies located beyond 
the virial radius of a cluster. This method can be applied even to t hose clusters w hich 


h ave yet to be completely relaxed without having spherical symmetry ( Corbett Moran et al. 


2 0141 ). However, large uncertainties inherently involved in the measurements of the peculiar 
velocity field remain as a technical bottleneck that stymies the accurate measurement of 
the dynamic mass of a cluster from the peculiar velocities of the neighbor galaxies (s ee also 
Corbett Moran et ah 2014 : Gronke et ah 2014b : Hellwing et al. 2014 : Zu et al. 2014 ) 


Falco et ah 1 20141) devised a novel method to measure the dynamic mass of a cluster by 
using the universal radial velocity profile of the filament galaxies around the cluster. The 
sole condition for the practical success of their method is that the neighbor galaxies should 
be distributed along a narrow filament. Since it has been theoretically and observationally 
shown that the for mation of galaxy clusters preferentially occur at the junction of the cosmic 
filaments (e.g., see lBond et al.lll996t iDietrich et al.ll20121) . the majority of the galaxy clusters 
sh ould meet this condition, being eligible for the application of this method. The method 
of lFalco et al. (20141) to measure the virial mass of a cluster from the radial velocity profile 
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of the filament galaxies is applicable even to those clusters which have yet to complete its 
relaxation process. In this respect their method fits especially well for the Virgo cluster 
that still undergoes merging events. Furthermore, the measurements of the redshifts and 
equatorial coordinates of the galaxies in the vicinity of the Virgo cluster were done with high 
accuracy thanks to their proximity, wh ich wi ll help reduc e uncertainties in the dynamical 
mass measurement with the method of Falco e t al. (120141 ) for the Virgo cluster. The goal 
of this paper is to test this new method against th e V irgo cl uster a round which the narrow 
filamentary structures were recently identified by Kim ct ah (120151 ). 


The upcoming sections will present the following. A review of the method of F a lco et al. 


1 20141 ) in Section [2] A reconstruction of the radial velocity profile of the filament galax¬ 
ies around the Virgo cluster and a new dynamical mass estimate for the Virgo clus ter by 


comparing the reconstructed profile with the numerical formula from Falco et al . (20141) 
in Section [3] Discussions on the final results and on the future works for the improve¬ 
ments in Section [4] Throughout this paper, we assume a GR+ACDM cosmolo gy wit h 


Q m = 0.26, Da = 0.74, h = 0.73 and zero curvature, to be consistent with F alco et al. 

(koi4lh 


2. METHODOLOGY : A BRIEF REVIEW 


The galaxies located in the vicinity of a cluster must feel the gravity of the cluster even 
in the case that they are not the satellites bound to the cluster, and thus their motions 
relative to the cluster should deviate from the pure Hubble flows. To quantify the degree 
of the deviati on o f th e g a laxy motions from the Hubble flows and its dependence on the 
cluster mass, Falco et_al.J (120141) determined the mean radial velocity profile of the galaxies 
in the vicinity of the clusters from a high-resolution N-body simulation for the GR+ACDM 
cosmology. Their numerical experiment has discovered that the mean radial velocity profile 
of the neighbor galaxies located in the zone where the effect of the Hubble expansion exceeds 
that of the gravitational attraction of a cluster is well approximated by the following universal 
formula: 

/ r \ 7i v 

v r (d-,z,M v ) = H(z)d- 0.8K V [jJ , (1) 

where v r is the mean radial velocity at (comoving) separation distance d from the cluster 
center, M v , r v and V v are the virial mass, radius and velocity of the cluster, respectively, and 
H(z) is the Hubble parameter at redshift z. The first term in the right-hand side represents 
a pure Hubble flow from the cluster center while the second term corresponds to the mean 
peculiar velocity induced by the gravitational attraction of the cluster. 


According to lFalco et al .l f 2 0141) . the radial velocity profile is universal in the respect that 
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the power-law index, n v , in Equation (JT|) has a constant value of 0.42, being independent of 
M v and z. They also showed that the zone where the radial velocities of the galaxies are well 
approximated by Equation dTj) with n v = 0.42 corresponds to the range of 3r v < d < 8r v . 
Adopting the conventional definition that r v is the radius of a sphere within whi ch the mass 


densit y reaches A c times the critical density of the Universe p c = 3H 2 /(8tt G), IFalco et al. 


12 01411 used the following simple relations to calculate the virial mass M v and the virial 
velocity V v for Equation (jT]): 


47T 

M v = —Ap c r^, 


U 2 = 


GM V 


( 2 ) 


The value of A c was set at 93.7 i n accordance w ith the formula of A c = 187T 2 + 82[O m (z) 
1] — 39[n m (js;) — l] 2 given by Brvan fc Norman ( 1998 k 


Falco et ah (120141 ) suggested that the virial mass of a massive cluster M v be estimated 


by comparing the observed radial velocity profile of its neighbor galaxies to Equation (JT]). 
However, since the real values of v r and d are not directly measurable, it is difficult to deter¬ 
mine the radial velocity profiles of the galaxies from observational data. They claimed that 
this difficulty should be overcome by using the anisotropic spatial distribution of the galaxies 
belonging to a one-dimensional filament in the vicinity of a cluster. Three dimensional sep¬ 
aration distance d and the radial velocity v r of the filament galaxies can be estimated from 
two direct observables, the projected distance R in the plane of the sky and the linc-of-sight 
velocity ui os , as R = sin (3r and v r = v\ os / cos/3, respectively, where f3 is the angle at which 
the filament is inclined to the line of sight direction to a cluster (see Figure [1]). Finally, 
Equation (pQ) was rewritten in terms of the observables as 


v los (R,(3,M v ) = cos/3 


H{z) 


R 

sin /3 


0.8K 


R 


sin /3 r v 


(3) 


Since the galaxy cl uster s usually form at the junction of the cosmic filaments, this new 
scheme of Falco et al. (120141 1 seems quite plausible. However, the validity of Equation (l3l) 
depends strongly on the geometrical shape of the filament. The more straight and narrow 
a filament is, the better Equation (JTJ) is approximated by Equation (J3J) . Figure [1] illustrates 
a three dimensional configuration of a filamentary structure around a cluster inclined at an 
angle (3 to the line of sight direction to the cluster, defining the projected d i stance R between 
a filament galaxy and the cluster center in the plane of sky. Falco et ah ( 120141 ) performed 
a numerical test of Equation (J3J) by applying it to three cluster-size halos around which 
filamentary structures were found. The estimated values of M v turned out to agree very well 
with the true virial masses for all of the three cases, which numerically verified the validity 
of Equations (J3J). 
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Falco et ah (2 014 ) also performed an observational test of their algorithm by applying 
it to the sheet structures identified around the Coma cluster. Despite their claim that the 
dynamic mass of the Coma cluster determined by their method turned out to be consistent 
with the previous estimates, their method, strictly speaking, is valid only for the galaxies 
distributed along a one-dimensional filamentary structure that interacts gravitationally with 
the cluster. In fact, it is obvious from Equation (j3J) that large uncertainties should contami¬ 
nate the measurement of the cluster mass if the radial velocity profile is reconstructed from 
the sheet galaxies. 


3. ESTIMATING THE DYNAMIC MASS OF THE VIRGO CLUSTER 


3.1. The Zero-th Order Approximation 


Kim ct ah (120151 ) developed an efficient technique to identify a filamentary structure 
whose galaxies are under the gravitational influence of a neighbor cluster. The merit of 
their filament-finding technique is that it requires information only on observable redshifts 
and two dimensional positions of the galaxies projected onto the plane of sky. Ap plying this 
fil ame nt-finding technique to the NASA-Sloan-Atlas catalog of the local galaxie^j, 


Kim ct ah 


(120151 ) have identified two narrow filamentary structures to the east and to the west of the 
Virgo cluster (say, the Filament A and B, respectively). Both of the Filament A and B were 
found to consist mainly of the dwarf galaxies with absolute r-band magnitudes M r > — 20. 
For a detailed descriptio n about how the Filament A and B were identified from the NASA- 
Sloan-Atlas catalog, see Kim ct ah (120151 ). 


Figure [2] plots the right asensions (R.A.) and declinations (decl.) of the galaxies be¬ 
longing to the filament A and B, as filled blue and red dots, respectively. The certain and 
possible mem ber g alaxies of the Virgo cluster from the Extended Virgo Cluster Catalog 
(EVCC) (IKim ct ahj 2014) are shown as the black and gray dots, respectively, inside an open 
rectangular box. The gray dots outside the box represent the field galaxies that do not belong 
to the EVCC nor to the two filaments but in the redshift range of 500 < cz / [km s _1 ] < 2500. 
Note that the Filament A appears to be more straight than the Filament B in the equatorial 
plane. 

Showing that the mean radial velocities of the ga laxies belonging to each filament are 
similar to the system velocity of the Virgo cluster, Kim ct ah (120151) have confirmed that the 
filament galaxies are in gravitational interaction with the Virgo cluster. Figure [H] plots the 


^iven in http://www.nsatlas.org/data 
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redshifts of the galaxies belonging to the Filament A and B as a function of their projected 
distances in unit of degree from the M87 (i.e., the Virgo center) as filled blue and red dots, 
respectively. As can be seen, the redshift variation of the filament galaxies A cz at the 
constant R are less than 1000 km/s. 


Denoting the equatorial positions (i.e., R.A. and decl.) of a filament galaxy and the 
Virgo center as (a g , S g ) and (a c i, 5 c i), respectively, we determine the projected distance R 
between each filament galaxy from the Virgo center in the plane of sky (see Figured]) as R — 
cz c i (cos 9 g cos 9 c i + sin^gsin^ci cosA0 g ) where 6 g = 7t/2 — S g , 6 C \ = 7t/ 2 — «5 c i, A0 g = a g — a c \, 
and z g and z& are the redshifts of the filament galaxy and the Virgo center, respectively. 
The line-of-sight velocity magnitude of each filament galaxy at a projected distance, R , 
from the Virgo center is approximated as c| z g (R) — z c \ |. 


Equation (J3J) is valid only for those filament galaxies whose separation distances from 
the clu ster center ar e larger than three times the virial radius of the Virgo cluster according 


to F alco et ah 020141 ). Taking the typical cluster size 2 h _1 Mpc as a conservative upper limit 
on the virial radius of the Virgo cluster, we select only those hlament galaxies which satisfy 
the condition of R > 6 /i _1 Mpc. A total of 130 and 96 gal axies ar e selec ted in the Filament A 
and B, respectively, to each of which the method of Falco et ah ( 20141 ) is applied as follows. 


Setting the power-law index n v at the universal value of 0.42, we compare the observed 
values of of the hlament galaxies with Equation (J3]). Assuming a hat prior, we search 
for the best-ht values of M v and (3 which maximize the following posterior density function 
p(M v , /31 data): 


p(M Vl fj\v° os , Ri ) oc exp s 


Nn 


i=l 


[v° os (Ri) - vio S {Ri] M v , 
2aUN g - 2 ) 


(4) 


where N g is the number of the selected hlament galaxies, Ri is the projected distance of 
the f-th hlament galaxy, cr v represents errors involved in the measurement of vtW from 
the redshift space. Given that cr v include all systematic errors that must be existent but 
unknown, we set cr v at unity for all of the selected hlament galaxies. Here, the poste¬ 
rior probability density function p(M v , (3\v° os , Ri) is normalized to satisfy the condition of 
f dM v f d/3p(M v , /3|uj° 0 g, R t ) = 1. 


The 68%, 95% and 99% confidence level contours of logM v and j3 obtained from the 
above fitting procedure for the case of the Filament A (B) are shown as solid, dashed and dot- 
dashed lines, respectively, in the left (right) panel of Figure U which reveals the existence of 
a strong degeneracy between M v and j3. To break this degeneracy, we would like to put some 
prior on the range of j3 by making the zeroth-order approximation to the three dimensional 
positions of the hlament galaxies relative to the Virgo cluster. 
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In the zero-tli order approximation of d ~ d° = cz g under the assumption of no pe¬ 
culiar velocity, the three dimensional separation vector of each selected filament galaxy, 
d° = (d\l, <7°, ^ 3 )) can be expressed as d® = c z g sin 8 g cos a g — cz c \ sin 5 c i cos a c \, d ° = 
c z g sin 5 g sin a g — c z c i sin d c i sin q c] , d 3 = cz g cos 5 g — cz c iCOs/) c i. The zeroth-order approxi¬ 
mation to the angle, /?, between d and the line of sight direction to the Virgo cluster, say h, 
is now calculated as, /3 ~ { 3 ° = cos _ 1 (d° • h) with d° = d°/|d°| where the ensemble average 
is taken over the selected filament galaxies. 


The Filament A and B yield /3° = 36.8° ± 3.9°and /3° = 36.9° ± 5.6°, respectively. It 
is interesting to note that this result is in line with the previous observational evidences for 
the alignm ent ten de ncy b etween the major axis of the Virgo cluster and the line-of-sight 
direction flMei et al.l 120071) . It is also consistent with the recent detection that the Virgo 
satellites tend to fall into the Virgo along the lie of s ight directions aligned with the minor 
principal axis of the local velocity shear held flLee et al.l 2014). Taking these zero-th order 
values as a prior, we break the degeneracy between M v and /3 shown in Figure [4] and finally 
determine the dynamical mass of the Virgo cluster: M v ps (0.841q^) x 10 15 h~ 1 M & for the 

(3. 215 ) 3 ) x 10 15 h~ l M Q for the case of the Filament B. 


case of the Filament A and M v ps 


3.2. Radial Motions of the Filament Galaxies around the Virgo Cluster 

It should be worth comparing our result with th e prev ious measurements of the Virgo 
mass based on various methods. For instance, Urban et al. (1201111 converted the X-ray 
temperature of the ho t gas i n the Virgo cluster to its mass via the mass-temperature scaling 
relation obtained by Arnaud et a l. ( 20051 1 and found M 200 /M Q fa 1.4 x 10 14 . Here M 2 00 
represents the mass enclosed by a spherical radi us r 2m at w hich the mass density becomes 
p(r 200 ) = 200 p cr it where p cr u = 3H 2 /8tiG. [Karachentsev et ah (2014j) studied the infall 
motions of seven nearby galaxies located along the line-of-sight direction to the Virgo cluster 
whose di stances were measured by means of the Tip of the Red Giant Branch (TRGB) 
method ( Baade 1944 : Salaris fe Cassisi 1997 1 and found M 2 oo/M Q = (7.0 ± 0.4) x 10 14 . 
It was also claimed that the lensing mass of the Virgo cluster shoul d be in the range of 
log (M 20 o/M q ) = 13.87l)) ) ;g5 by employing the technique developed by Kubo et al . ( 2009 1 to 
estimate the lensing mass of the galaxy clusters in the local universe (D. Nelson 2014, private 
communicatio .ft. 


Noting that the definition of M 2 oo is different from that of the virial mass M v considered 
here, we convert the values of M 2 qq/Mq reported in the previous works into M v / (hr 1 M 0 ) 


2 see also the presentation given in https://www.cfa.harvard.edu/ dnelson/doc/dnelson.fermilab.talk.pdf. 




























for a fair comparison under the assumption that the mass d ensity profile of the Virgo cluster 
is well approximated by the NFW formula (INavarro et al.l 119971 ) and list them in Table [T] 
As can be read, our best-fit value M v obtained from the Filament A is consistent with the 
previous dynamic mass estimated by the TRGB-distance method but substantially higher 
than the other two estimates. Meanwhile, for the case of the Filament B, the best-fit M v 
is significantly off not only from the previous estimates bud from the best-fi t value obtained 
for the case of the Filament A. Given that the method of [Falco et ah (12014 ) is strictly valid 
only for the galaxies distributed along thin straight filament, we suspect that the bent shape 
of the Filament B should be responsible for the significant discrepancy on M v between the 
two filaments and suggest that the Filament B should not be eligible for the application of 
the method of Falco et al. ( 2014 ). 


Before excluding the Filament B from our analysis, it may be worth investigating how 
different the reconstructed velocity pro file of th e Virg o cluster is from the expected profile, 
i.e., the universal formula given by [Falco et ah (12014 ) when the mass of the Virgo cluster is 
fixed at some constant value M v . Treating n w in Equation (j3J) as a free parameter, we search 
for the best-fit values of n v and /? that maximize the following posterior distribution: 


p(n v ,/3\v£ s , Ri) oc p(data|n v ,/3)p(/3) 


oc exp 


f Kb (Ri) - Vlos {Ri\n v , /3)] 2 ' 

1 . 

1 

CM 

1 

1 _ 

2 (N, - 2) J 

> exp 

2o^ 0 


(5) 


Here, instead of setting (3 at the fixed value (3°, we assume a Gaussian prior on (3, setting 
its mean and standard deviation at (3° and a^o, respectively, both of which are determined 
in Section 13.11 


Figure 0 shows the 68 %, 95% and 99% confidence level contours of f3 and n v obtained 
from the Filament A as solid, dashed and dot-dashed lines, respectively, for the four different 
cases of the dynamical mass of the Virgo cluster: M v = 5 x 10 13 , 10 14 , 5 x 10 14 , 10 15 /i _ 1 M 0 . 
As can be seen, for the cases of M v < 10 14 h _ 1 M 0 that is consistent with the previous 
estimate of the lensing mass (see Table 1), the original value of n v = 0.42 is located beyond 
the 95% confidence level contours. For the case of M v > 5 x 10 14 h~ 1 M & that is consistent 
with the dynamical mass estimate obtained by the TRGB method, n v = 0.42 is included 
within the 68 % confidence level contour. H ence, one can conc lude that the numerical formula 
of the radial velocity profile suggested by Falco et al. (120141 ) matches fairly well with the 
reconstructed prohle from observation for the case of the Filament A around the Virgo cluster, 
yielding the best-fit value of M v consistent with the previous dynamical mass measurement 
for the Virgo. Figure [ 6 ] plots the same as Figure[5]but from the Filament B. As can be seen, for 
all of the four cases of M v , the original value of 0.42 stays outside the 95% contours, indicating 
that Equation ([3]) fails to describe the radial motions of the galaxies in the Filament B. 























9 


We also show the radial velocity profile, v r (d), reconstructed from the Filament A as 
black solid line in Figure [7] To plot v r , the Virgo mass M v is set at 10 14 h~ 1 M@ (close to the 
lensing mass), the slope n v and the angle (3 are set at the best-fit values determined by the 
above fitting procedure. The associated la and 2a uncertainties around the reconstructed 
profile are also shown as dark and light grey regions, respectively, and the red dashed line 


r epres ents the expected profile with the fixed value of n v = 0.42 claimed by iFalco et ah 


(12014T) . As can be seen, the expected profile is only marginally consistent with the recon¬ 
structed one if M v = 10 14 /i" 1 M 0 . Figure [H] plots the same but for the case of the Filament 
B. It is clear that the difference between the numerical profile and the reconstructed one 
becomes much larger for the case of the Filament B. 


4. DISCUSSION AND CONCLUSION 


We have measured the dynamic mass of the Virgo cluster, M v , by applying the method 
of F alco e t al. (1201 4) to two narrow filaments (the Filament A and B) consisting mainly of 
dwarf galaxies located to the east and to the west of the Virgo cluster in the equatorial 
reference frame, which were r ecen tly identi fied b y Kim et. ah ( 20151) from the NASA-Sloan- 
Atlas catalog. The method of lFalco et al.l 12 01411 that can be applied even to the unrelaxed 
clusters like the Virgo is based on the numerical findings that the radial velocity profile of 
the galaxies at distances larger than three times the virial radius of a neighbor cluster has 
a universal shape and that it can be readily reconstructed from direct observables as far as 
the galaxies are distributed along a filamentary structure. 

To break a strong degeneracy found between the value of M v and the angle /?, we 
have constrained the ranges of /3 by making a zero-th order approximation in which the 
redshift-space positions of the filament galaxies are assu med to equal their real-space ones. 
Our estimates based on the method of Fa lc o et al. ( 2014 ) has yielded M v = (0.84lgl5i) x 


10 15 h 1 M q and M v = (3.2i^g) x 10 15 h 4 M 0 from the Filament A and B, respectively. 

Given that the Filament B appears not so straight as the Filament A in the equatorial 
plane, the substantial difference between the two results on the value of M v from the Fila¬ 
ment A and B h as been suspected as an indication that the Filament B is not eligible for 
the algorithm of lFalco et al . (2014) wh ose success depen ds sensitively on the degree of the 
filament straightness. The method of Falco et al. (20141) itself also suffers from a couple of 
unjustified assumptions. First of all, it assumes the filaments to be extremely narrow and 
straight, which is an obvious over-simplication of the true shapes of the filaments. Second 
of all, in their original study, the universal radial velocity profile was obtained by taking an 
average over all the particles and halos around the clusters but not only over the filament 
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halos. If only a subsample composed of the filament galaxies is used to reconstruct the radial 
velocity, then it would not necessarily agree with the average universal profile. Besides, the 
zero-th order approximation that we have made for f3 in our estimation of the Virgo mass 
with the method of IFalco et a l. ( 20141 ) is an additional simplification of the reality. If a 
broader range of f3 were considered, then the mismatches between the two masses would be 
likely reduced. 

Our estimate of M v from the Filament A has turned out to be consi stent w ith th e prev i¬ 


ous dy namical mass estimates based on the TRGB-distance method by [Karachentsev et al. 


12 0141 1 but three to four times higher than the other mass estimates based on the weak gravi¬ 
tational len sing effect fD . Nelson 2014 in private communication) and the X-ray temperature 
of hot gas flUrban et al. 201 111. The usual suspect for this disagreement between the mass 
estimates for the Virgo cluster is the contamination of the measurements due to the presence 
of systemati cs associated w ith such over-simplified assumptions as the spherical NFW den¬ 
sity profile (Navarro et al. 19971 ) in the lensing mass estimates, the hydrostatic equilibrium 
condition for the X-ray temperature measurement, and etc. 

Another more speculative explanation for the difference between the lensing and the 
dynamic mass measurements for the Virgo cluster is the failure of GR. We have shown that if 
the dynamic mass of the Virgo cluster is assumed to be identical to the lensing mass, then the 
reconstructed radial velocity profile from observational data appears to less rapidly decrease 
with the distance than the universal radial profile from the simulation for the GR+ACDM 
cosmology. The less rapid decrease of the radial velocity profile with distance indicates higher 
peculiar velocities of the filament galaxies at large distances than predicted for the case of 
GR. In fact, several numerical works already explored how much the peculiar velocities of 
galaxies around massive clusters would be enh a nced in the presence of modified grav ity (MG) 
( Lam et al. 2012 ; Corbett Moran et al. 2014 ; Gronke et al. 2014b ; Zu et al. 2014). Before 
exploring this speculative idea, however, it will be necessary to study comprehensively how 
the systematics produced by the over-simplified assumptions mentioned in the above would 
affect the reconstruction of the radial velocity profile of the filament galaxies in the vicinity 
of the clusters. Our future work is in this direction. 
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Fig. 1.— Illustration of a three dimensional configuration of a filamentary structure inclined 
at an angle /? to the Virgo cluster. The projected distance R of a filament galaxy from the 
Virgo center in the plane of sky perpendicular to the line of sight direction toward the Virgo 
can be measured from information on their equatorial coordinates. 
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Fig. 2.— Spatial distribution of the Filament A and B around the Virgo cluster in the 
equatorial coordinate system. The black dots correspond to the galaxies of the Virgo Cluster 
while the filled blue (red) dots correspond to the galaxies belonging to the Filament A (B). 
The grey dots represent the galaxies that are not associated with the Virgo cluster nor 
filaments.The fitted lines of the galaxies in the Filament A and B are also denoted by the 
solid curves. 
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Fig. 3.— Distances of the Filament galaxies from the center of the Virgo cluster versus their 
recession velocities. The symbols are the same as in Figured 
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Fig. 4.— 68% (solid), 95% (dashed) and 99% (dot-dashed) confidence level contours of logM v 
and /3 for a ACDM cosmology by comparing Equation (J3]) to the observed radial velocity 
profile of the galaxies in the Filament A and B in the left and right panels, respectively. In 
each panel the cross marks the best-fit (log M v ,/3). 
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Fig. 5.— 68% (solid), 95% (dashed) and 99% (dot-dashed) confidence level contours of n v 
and f3 from the observed radial velocity profiles of the galaxies in the Virgo Filament A. The 
cross marks the best-fit (n v ,(3). 
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Fig. 6.— Same as Figure [5] but for the case of the Filament B. 
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Fig. 7.— Velocity profile of the Filament A galaxies around the Virgo. The black solid line 
corresponds to the best-fit velocity profile inferred from the observation data while the red 
solid line is the expected velocity profile of a galaxy cluster with mass of 10 14 h~ 1 M & in the 
ACDM cosmology. The regions in dark (light) gray colors represent the la (2 a) scatters in 
the determination of the best-fit velocity profile 
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Fig. 8.— Same as Figure [T] but for the case of the Filament B galaxies. 
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Table 1. Observable, virial mass estimate of the Virgo cluster in unit of 10 14 h 1 M 0 with 
h = 0.73 and literature that reports the estimate. 


indicator 

M v 

literature 


(10 14 H-'Mq) 


X-ray temperature 

1.50 

Urban et al. (2011) 

weak lensing effect 

o.8ii5;g 

Nelson (2014) a 

TRGB distance 

7.73 ±0.43 

Karachentsev et al. (2014) 


a https://www.cfa.harvard.edu/ dnelson/doc/dnelson.fermilab.talk.pdf 






